Stability of the Einstein static universe in modified Gauss-Bonnet gravity 
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We analyze the stability of the Einstein static universe by considering linear homogeneous pertur- 
bations in the context of f(G) modified Gauss-Bonnet theories of gravity. By considering a generic 
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equation of state parameter w = p/p and the second derivative f"(G) of the Gauss-Bonnet term. 
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I. INTRODUCTION 

Cosmology is said to be thriving in a golden age, where 
a central theme is the perplexing fact that the Universe is 
undergoing an accelerating expansion The latter, one 
of the most important and challenging current problems 
in cosmology, represents a new imbalance in the govern- 
ing gravitational equations. Historically, physics has ad- 
dressed such imbalances by either identifying sources that 
were previously unaccounted for, or by altering the gov- 
erning equations. The cause of this acceleration still re- 
mains an open and tantalizing question. Although the in- 
troduction of a cosmological constant into the field equa- 
tions seems to be the simplest theoretical approach to 
generate a phase of accelerated expansion, several alter- 
native candidates have been proposed in the literature, 
ranging from dynamical dark energy models to modified 
theories of gravity. Amongst the latter, models general- 
izing the Einstcin-Hilbert action have been proposed. 

The Einstein field equations of General Relativity were 
first derived from an action principle by Hilbert, by 
adopting a linear functional of the scalar curvature, R, 
in the gravitational Lagrangian density. However, there 
are no a priori reasons to restrict the gravitational La- 
grangian to this form, and indeed several generaliza- 
tions of the Einstein-Hilbert Lagrangian have been pro- 
posed, including "quadratic Lagrangians" , involving sec- 
ond order curvature invariants such as R 2 , i? ji4 ,i? A " / , 
R a ^R af3tlu , e a ^R a0lS R^i, C a0IMV C a ^», etc (see 
Ref. 3 for a recent review) . The physical motivations for 
these modifications of gravity were related to the possibil- 
ity of a more realistic representation of the gravitational 
fields near curvature singularities and to create some first 
order approximation for the quantum theory of gravita- 
tional fields. 

In considering alternative higher-order gravity theo- 
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ries, one is liable to be motivated in pursuing models 
consistent and inspired by several candidates of a fun- 
damental theory of quantum gravity. In this context, it 
may be possible that unusual gravity-matter couplings 
predicted by string/M-theory may become important at 
the recent low-curvature Universe. For instance, one may 
couple a scalar field not only with the curvature scalar, 
as in scalar-tensor theories, but also with higher order 
curvature invariants. Indeed, motivations from string/M- 
theory predict that scalar field couplings with the Gauss- 
Bonnet invariant Q are important in the appearance of 
non-singular early time cosmologies [H ■ It is also possible 
to apply these motivations to the late-time Universe in 
an effective Gauss-Bonnet dark energy model 0. 

In the context of f(Q) modified theories of gravity, 
we explore the stability of the Einstein static universe 
in this work. This is motivated by the possibility that 
the universe might have started out in an asymptoti- 
cally Einstein static state, in the inflationary universe 
context 0. On the other hand, the Einstein cosmos has 
always been of great interest in various gravitational the- 
ories. In general relativity for instance, generalizations 
with non-constant pressure have been analyzed in [3] . In 
brane world models, the Einstein static universe was in- 
vestigated in [8| , while its generalization within Einstein- 
Cartan theory can be found in Hfl, and in loop quantum 
cosmology, we refer the reader to [10| . 

In the context of f(R) modified theories of gravity, the 
stability of the Einstein static universe was also analyzed 
by considering homogeneous perturbations [ll|. By con- 
sidering specific forms of f(R), the stability regions of the 
solutions were parameterized by a linear equation of state 
parameter w — p/p. Contrary to classical general rela- 
tivity, it was found that in f(R) gravity a stable Einstein 
cosmos with a positive cosmological constant does indeed 
exist. Thus, in principle, modifications in f(R) gravity 
stabilize solutions which are unstable in general relativ- 
ity. Furthermore, in [12| it was found that only one class 
of f(R) theories admits an Einstein static model, and 
that this class is neutrally stable with respect to vector 
and tensor perturbations for all equations of state on all 
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scales. These results are apparently contradictory with 
those of Ref. [ll|. However, in a recent work, homo- 
geneous and inhomogeneous scalar perturbations in the 
Einstein static solutions were analyzed (l3j , consequently 
reconciling both of the above works. 

This paper is outlined in the following manner: In Sec- 
tion [Til we briefly review modified Gauss-Bonnet gravity, 
for self-completeness and self-consistency, and present 
the respective field equations. In Section IIII1 we con- 
sider a generic form of f(Q), and analyze the stability of 
the solutions, by considering homogeneous perturbations 
around the Einstein static universe, and deduce the re- 
spective stability regions. The latter are given in terms 
of the linear equation of state parameter w — p/ p and 
the unperturbed energy density po- Finally, in Section 
IIVI we present our conclusions. Throughout this work, 
we consider the following units c = G = 1 . 



II. MODIFIED GAUSS-BONNET GRAVITY 
AND FIELD EQUATIONS 

An interesting alternative gravitational theory is modi- 
fied Gauss-Bonnet gravity, which is given by the following 
action 



S = d 4 x 



+ S M (g f iu,ip), (i) 



where k = 8n. f(G) is an arbitrary function of the Gauss- 
Bonnet invariant, which is in turn given by 



g = r 2 - m^R^ + Ru Va pR^ ap . 



(2) 



Sm(9hv,4>) is the matter action, defined as Sm = 
J d x-\/—g £ m (g M „, V0> where C m is the matter La- 



grangian density, in which matter is minimally coupled 
to the metric g pv and ip collectively denotes the matter 
fields. The matter stress-energy tensor, , is defined 



as 



Thus, using the diffeomor- 



phism invariance of Suig^v, ip) yields the covariant con- 



servation of the stress-energy tensor, V^T^ 



(m) _ 



0. 



Let us for the moment assume the function f(G) to be 
expanded around the origin so that we may write 



f(G) = /(o) + f(o)g 



no) 

2! 



g 2 + o(g 3 ). (3) 



If we insert this representation back into the Lagrangian, 
we note that /(0) acts like an effective cosmological con- 
stant. The term /'(0) would not be present in the field 
equations because this term is linear in the Gauss-Bonnet 
invariant and linear topological terms do not contribute 
to the equations of motion. The first non-trivial term to 
be expected is f"(0). Moreover, as we are interested in 
linear perturbation theory, we expect only the term /"(0) 
to characterize deviations from general relativity. Hence, 
we anticipate that the general relativistic limit takes the 
form /"(0) -► 0. 

Modified Gauss-Bonnet gravity has been extensively 
analyzed in the literature, and rather than review all 
of its intricate details here, we refer the reader to 
Refs. HI EE El- Varying the action Eq. dJ) with re- 
spect to g 1 *", one obtains the gravitational field equa- 
tions in the following form G» v = R» v - \g^ v R = kT^, 
where the effective stress-energy tensor is defined as 
^eff = r ^(m)~^ r ^GB- The effective Gauss-Bonnet curvature 
stress-energy term is given by 



J GB - 2 



<r7(S) - 2f(g)RR^ + 4f XQ)R» p R v p - 2f'(g)R^ T R% aT ~ ^'{g)R^ av R pa 



- 2g^[w 2 f(g)]R - iiVpV^f'igw - , ... 

+ Ag^[V p V a f'{g)}R p ' y - 4[V p V CT /'(S)]i? WWCT 



4[V P V7'(£)]-R W 



4[V7'(6)]fl"'' 



2{V^V y f(g)]R 



(4) 



where the definition /'(£?) 
notational simplicity. 



df(g)/dg is considered for III. 



THE EINSTEIN STATIC UNIVERSE IN f(Q) 
GRAVITY AND PERTURBATIONS 



Metric and field equations 



Consider the metric given by 



ds 2 = -dt 2 



dr 2 
1 — r 5 



+ r 2 (d9 2 +sm 2 9d(/) 2 ) 



(5) 
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The Gauss-Bonnet invariant in provided by 



Q = 24 



1 + a! 1 a" 



(6) 



where the prime, in this context, denotes differentiation 
with respect to cosmological time t. 

For the Einstein static universe, a = ao = const, the 
Ricci scalar becomes R = 6/0$, and the Gauss-Bonnet 
invariant is trivially given by Q(ao) = 0. The field equa- 
tions in this case take the following form 



3 ( /(0) 



a, 



where po and po are the unperturbed energy density 
and isotropic pressure, respectively. It should be noted 
that the modified Gauss-Bonnet term /(0) acts like an 
effective cosmological constant to the background field 
equations, as noted above. 



The first equation relates the perturbations in the scale 
factor to the density perturbations. Next, we can elimi- 
nate 5p from the second equation using the first and we 
arrive at the following fourth order perturbation equa- 
tions for the perturbed scale factor 

2 f"(0) 
-*(1 + 3w)5a - 25a" - Wn—^-Sa^ = . (13) 

By virtue of the background equations we have 



—s = + W)p[ 







(14) 



and therefore we finally arrive at the following perturba- 
tion equation 

24 K 3 po(l + H 2 /"(0)<W 4) W 

+ 25a(t) - np (l + w)(l + 3w)Sa(t) = . (15) 



B. Linear homogeneous perturbations 



C. Stability regions 



In what follows, we analyze the stability against linear 
homogeneous perturbations around the Einstein static 
universe given in Eqs. ([7]). Thus, we introduce pertur- 
bations in the energy density and the metric scale factor 
which only depend on time 

p(t) = po (1 + 8p{t)) , a{t) = a (1 + Sa(t)) . (8) 

Subsequently, we consider a linear equation of state, 
p(t) = wp(t), linearize the perturbed field equations and 
analyze the dynamics of the solutions. 

Despite the Gauss-Bonnet invariant being trivial at the 
background level, the perturbations have an effect. From 
Eq. © it follows that 



% a 



(9) 



in linear order. Therefore, using f(Q + 5Q) = f(Q) + 
f'{Q)8Q plus higher order terms, we find for linear per- 
turbations 



h 8n" 

/(0) = /(O) + 24/'(O)4^-, 
% a 



(10) 



where we assumed that / is an analytic function. 

We now insert the perturbations |5]) into the full field 
equations and linearize. This provides the following two 
equations 

-^8a + Kp o 5p = (11) 
a o 

2 f"(Q) 
—5a - 25a" - K Po w5p - 96k— ^fo (4) = (12) 



Using the standard ansatz 5a(t) — Cexp(wi), where 
C and to are constants, we find that Eq. (fT"5|) provides 
as solutions the following four frequencies of the small 
perturbations 



-1 ± y/1 + 24« 4 ^(1 + w) 3 (l + 3w)/"(0) 
24k 3 / 9q( 1 + w) 2 /"(0) 



(16) 



As expected, in the limit /"(0) — > the frequencies u)\ 
lead to the general relativistic result 
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npo 



(l + tfl)(l + 3tD) 



(17) 



while the frequencies to^_ become formally infinite in this 
limit, which is simply an artifact of the reduction of 
fourth-order to second order gravity. 

The perturbations 5a are stable if the frequencies uj 
are purely complex, which implies that the stability is 
equivalent to the two conditions u>\ < 0. In order to 
simplify the notation, let us introduce the new parameter 
a = 24k 4 /5q/"(0), so that the perturbations are stable for 
the following three parameter regions. 

First, there is a region analogous to the classical gen- 
eral relativistic stability region given by 



i >w> 



1 



(l + w) 3 (l + 3w) 



>a>0. (18) 



The second region covers the range where the equation 
of state is larger due to the presence of the Gauss-Bonnet 
term, and is provided by 



It should be noted that these perturbation equations do 
not contain contributions of the form /'(0), as expected. 



1 

w>-~. 



> a > - 



1 



(l + w) 3 (l + 3w) 



(19) 
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Finally, there exits a stable region where the equation 
of state is less than minus one. We denote this the phan- 
tom region, and the stability regions are given by 



w < -1 , 



> a > 



(1 + w) 3 (l + 3w) 



(20) 



These three regions are depicted in Fig. [T] and it is 
evident that there exits stable modes for all equation of 
state parameters w if a is chosen appropriately. We em- 
phasize that this includes positive values of the equation 
of state. It is interesting to note that the negative values 
of a = 24/t 4 po/"(0), imply that /"(0) < 0. 




equation of state w 



FIG. 1: Regions of stability in the (w,a) parameter space 
for homogeneous perturbations of Einstein static universes. 
The left is the phantom region, the middle part is the general 
relativistic analogue region and the right part corresponds to 
a normal matter region. Note that there exist stable modes 
for all w, provided a is chosen appropriately. 



as an inflationary cosmology without a singularity |6j. 
The role of positive curvature, negligible at late times, 
is crucial in the early universe, as it allows these cos- 
mologies to inflate and later reheat to a hot big-bang 
epoch. An attractive feature of these cosmological mod- 
els is the absence of a singularity, of an 'initial time', of 
the horizon problem, and the quantum regime can even 
be avoided. Furthermore, the Einstein static universe 
was found to be neutrally stable against inhomogeneous 
linear vector and tensor perturbations, and against scalar 
density perturbations provided that the speed of sound 
satisfies > 1/5 [l7|. Further issues related to the sta- 
bility of the Einstein static universe may be found in 
Ref. [3. 

In this work we have analyzed linear homogeneous per- 
turbations around the Einstein static universe in the con- 
text of f(G) modified theories of gravity. In particu- 
lar, perturbations in the energy density and the metric 
scale factor were introduced, a linear equation of state, 
p(t) = wp(t), was considered, and finally the linearized 
perturbed field equations and the dynamics of the solu- 
tions were analyzed. It was shown that stable modes for 
all equation of state parameters w exist, if the parameter 
a is chosen appropriately. Thus, as in Refs. [IJ, [l]| 
our results show that perturbation theory of modified 
theories of gravity present a richer stability/instability 
structure that in general relativity. Finally, it is of inter- 
est to extend our results to inhomogeneous perturbations 
m the spirit of Ref. [lj], 

and to include the canonical 
scalar field case. Work along these lines is presently un- 
derway. 



IV. SUMMARY AND DISCUSSION 
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